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Abstract. After introducing the basic ingredients of Loop Quantum Cosmology, I will briefly 
discuss some of its phenomenological aspects. Those can give some useful insight about the full 
Loop Quantum Gravity theory and provide an answer to some long-standing questions in early 
universe cosmology. 



1. Introduction 

The variety of precise astro-physical and cosmological data, available at present, combined with 
the large number of high energy physics experiments, are expected to give us the necessary 
ingredients to test fundamental theories and understand the very early phases of the evolution 
of our universe. 

Cosmological inflation p] remains the most promising candidate to solve the shortcomings 
of the standard hot big bang model, while it offers a causal explanation for the primordial 
fluctuations with the correct Cosmic Microwave Background (CMB) features. Despite this 
success, inflation suffers from a number of drawbacks [2j. In particular, compatibility between 
theory and measurements often necessitates fine-tuning of the inflationary parameters, inflation 
remains still a paradigm in search of a model, while it must prove itself generic [21 IH El El El El E] • 

There is undoubtedly an additional list of fundamental cosmological questions, still lacking 
a satisfying answer. One does not know, for instance, how close to the big bang a smooth 
space-time can be considered as the correct framework. Quantum gravity, a full theory which is 
supposed to resolve the big bang singularity, is still missing, while it is not known whether a new 
boundary condition is needed at the big bang, or whether quantum dynamical equations remain 
well-behaved at singularities. It is nevertheless clear that a smooth space-time background 
cannot be assumed as the correct description close to the big bang. 

In a Hamiltonian formulation to a quantum theory, the absence of background metric 
indicates that Hamiltonian dynamics is generated by constraints. Physical states are solutions 
to quantum constraints implying that all physical laws are obtained from these solutions, while 
there is no external time according which evolution can be studied. One has to define a monotonic 
variable to play the role of emergent, or internal, time, and then build a framework within which 
the short-distance drawbacks of General Relativity near the big bang can be cured, maintaining 
however the agreement with General Relativity at large scales. 

Loop Quantum Gravity (LQG) |101 [TT] is a non-perturbative and background independent 
canonical quantisation of General Relativity in four space-time dimensions. Loop Quantum 
Cosmology (LQC) [12] is a cosmological mini-superspace model quantised with methods of 
LQG. The discreteness of spatial geometry and the simplicity of the setting allow for a complete 



study of dynamics. The difference between LQC and other approaches of quantum cosmology, 
is that the input is motivated by a full quantum gravity theory. The simplicity of the setting, 
combined with the discreteness of spatial geometry provided by LQG, render feasible the overall 
study of LQC dynamics. 

In what follows, I will briefly discuss some of the phenomenological consequences of LQC, 
a subject which gains constantly an increasing interest from the scientific community, due in 
particular to its successes. 

2. Elements of LQG/LQC 

Loop quantum cosmology is formulated in terms of SU(2) holonomies of the connection and 
triads. The canonically conjugated variables are defined by the densitised triad Ef, and an 
SU(2) valued connection A l a , where i refers to the Lie algebra index, and a is a spatial index 
with a, i = 1, 2, 3. The densitised triad gives information about spatial geometry (i.e., the three 
metric), while the connection gives information about curvature (spatial and extrinsic one). 
Let me explain this: Quantum gravity introduces a discreteness to space-time. To quantise 
quantum gravity, this discreteness manifests itself as quanta of space. Since the quanta are three- 
dimensional, they can be characterised by a triad of numbers. The connection between adjacent 
pairs of these quanta form a two-dimensional surface, through which the geometric connection is 
defined. To ensure that local rotations do not induce different geometries, this connection must 
be SO(3) invariant. By using connection-triad variables, arising from a canonical transformation 
of Arnowitt-Desner-Misner ( variables, one can make an analogy with gauge theories, which is 
particularly useful when dealing with quantisation issues. 

For any quantisation scheme based on a Hamiltonian framework (e.g., LQC) or an action 
principle (e.g., path integral) for a homogeneous and flat model, divergences appear which need 
to be regularised. To remove the divergences that occur in non-compact topologies, the spatial 
homogeneity and Hamiltonian are restricted to a finite fiducial cell, of scale with finite volume 
Vo = f d 3 xy^\°q\, where °q is the determinant of the fiducial background metric. 

In the case of a spatially flat background, derived from the Bianchi I model, the isotropic 
connection can be expressed in terms of the dynamical component of the connection c(t) as 



with u l a a basis of left-invariant one- forms ui l a = dx l . The densitised triad can be decomposed 
using the Bianchi I basis vector fields Xf = Sf as 



where °q stands for the determinant of the fiducial background metric q a b = oj^oju, and p(t) 
denotes the remaining dynamical quantity after symmetry reduction. 



where the absolute value is taken because the triad has an orientation. Since the basis vector 
fields are spatially constant in the spatially flat model, the connection component is 



K = c(t)< 
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c = sgn(p)7 
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N is the lapse function and 7 the Barbero-Immirzi parameter, a quantum ambiguity parameter 
approximately equal to 0.2375. 



The canonical variables c,p are related through the Poisson bracket 



{c,p} = ^V , (5) 

where Vq the volume of the elementary cell adapted to the fiducial triad and k = 8ttG. 

Defining the triad component p, determining the physical volume of the fiducial cell, and 
the connection component c, determining the rate of change of the physical edge length of the 
fiducial cell, as 

p = V Q 2/3 p , c = V 1/3 c, (6) 

respectively, we obtain 

{c, P } = ^, (7) 

independent of the volume Vq of the fiducial cell. 

To quantise the theory one faces the usual difficulty of quantum cosmology. Namely, the 
metric itself has to be considered as a physical field, thus it must be quantised; it is not a 
fixed background. Any standard quantisation scheme fails: Fock quantisation fails, while even 
free Hamiltonians are non-quadratic in metric dependence, so there is no simple perturbation 
analysis. The way out is to use gauge theory variables to define holonomies of the connection 
along a given edge 

h e (A) = Pex P y d S 7 M (s)4(7(s)M > (8) 

where V indicates a path ordering of the exponential, 7^ is a vector tangent to the edge and 
Ti = —i(Ti/2, with <7j the Pauli spin matrices, and fluxes of a triad along an S surface 

E(S, f)= I e abc E ci fidx a dx b , (9) 
Js 

with fi an SU(2) valued test function. Certainly, at the level of LQC these variables may seem 
ad-hoc and unnatural, however their motivation follows naturally from the full LQG theory. 
Thus, the basic configuration variables in LQC are holonomies of the connection 

h^\A) = cos (ff )l + 2sin(^)r 4 , (10) 

along a line segment ^0 ef and the flux of the triad 

F s {E,f)<xp, 



where the basic momentum variable is the triad component p, 1 is the identity 2x2 matrix and 
Ti = —iai/2 is a basis in the Lie algebra SU(2) satisfying the relation 

T iTj = (l/2)e ijk T k - (l/4)(% . 

In the classical theory, curvature can be expressed as a limit of the holonomies around a loop 
as the area enclosed by the loop shrinks to zero. In quantum geometry however, the loop cannot 
be continuously shrunk to zero area and the eigenvalues of the area operator are discrete. Thus, 
there is a smallest non-zero eigenvalue, the area gap A [13]. As a result, the Wheeler-de Witt 
differential equation gets replaced by a difference equation whose step size is controlled by A. 

Let me start with the old quantisation procedure. Along the lines of LQG, one considers 
e ifj.oc/2 (with ^ Q an arbitrary real number) and p, as the elementary classical variables, which 



have well-defined analogues. Using the Dirac bra-ket notation and setting e l ^° C//2 = (c|/i), the 
action of the operator p acting on the basis states reads 



PW = — Im> 



where \i (a real number) stands for the eigenstates of p, satisfying the orthonormality relation 



(M1IM2) = <Wm2" The action of the exp 



operator acting on basis states is 



exp 



«Mo 



\fi) = exp 



Mo 



|m) = Im + Mo) , 



(12) 



where /io is any real number. Thus, in the old quantisation, the operator e J/i ° C//2 acts as a simple 
shift operator. As in the full LQG theory, there is no operator corresponding to the connection, 
however the action of its holonomy is well defined. The action of the holonomies, h!f"°\ of the 
gravitational connection on the basis states is given by [H 



= (csl + 2 S nTi)|M> 



(13) 



where, 



cs|/x) = cos(/i c/2)|/x) 
m\n) = sin(fj,oc/2)\fj) 



[ Im + Mo) + Im~Mo) ]/2 , 
[ lM + Mo) - Im-Mo) ]/(2») • 



(14) 



The gravitational part of the Hamiltonian operator in terms of SU(2) holonomies and the 
triad component, in the irreducible J = 1/2 representation, reads |14j 
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sgn(p) . (15) 



The action of the self-adjoint Hamiltonian constraint operator, jH s = (C g + C g )/2 on the basis 
states \fj) is 



WglM) 



where 



{ + R(p + 4/i )] |m + 4/io) - 4i?(M) |m> + [^(m) + #(m - 4 Mo)] Im - 4mo> } , 

(16) 
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The dynamics are then determined by the Hamiltonian constraint 

(H g + H (t> m = , 



3/2 



(17) 



(18) 



where stands for the matter Hamiltonian. In the full LQG theory there is an infinite number 
of constraints, while in LQC there is only one integrated Hamiltonian constraint. Matter is 
introduced by just adding the actions of matter components to the gravitational action. We can 
finally obtain difference equations analogous to the differential Wheeler-de Witt equations. 



More precisely, the constraint equation on the physical wave- functions |^), which can be 
expanded using the basis states as |^) = ^V(0)Im) with summation over values of fi and 
where the dependence of the coefficients on represents the matter degrees of freedom, reads |15j 



+ 



V, 



+ 



+ 



^+4^(0) -4 



= - 4k2 ^ ^(m)^W ; (19) 



the matter Hamiltonian is assumed to act diagonally on the basis states with eigenvalues 
H^fi). Equation (fT9|) is the quantum evolution (in internal time fi) equation; there is no 
continuous variable (the scale factor in classical cosmology), but a label /i with discrete steps. 
The wave- function \I>„(0), depending on internal time [i and matter fields <f>, determines the 
dependence of matter fields on the evolution of the universe, with a massless scalar field playing 
the role of emergent time. Thus, in LQC the quantum evolution is governed by a second order 
difference equation, rather than the second order differential equation of the Wheeler-de Witt 
quantum cosmology. As the universe becomes large and enters the semi-classical regime, its 
evolution can be approximated by the differential Wheeler-de Witt equation. 



3. Lattice refinement 

The quantised holonomies were at first assumed to be shift operators with a fixed magnitude, 
leading to a quantised Hamiltonian constraint being a difference equation with a constant interval 
between points on the lattice. While these models can be used to study certain aspects of the 
quantum regime, they were found to lead to serious instabilities in the continuum, semi-classical 
limit |16l 117] . Considering the continuum limit of the Hamiltonian constraint operator, we 
have so far assumed that \P does not vary much on scales of the order of 4fj,Q (known as pre- 
classicality) , so that one can smoothly interpolate between the points on the discrete function 
^/i((p) and approximate them by the continuous function ^(/i, (f>). Under this assumption, the 
difference equation can be very well approximated by a differential equation for a continuous 
wave-function. However, the form of the wave-function indicates that the period of oscillations 
can decrease as the scale increases. Thus, the assumption of pre-classicality can break down at 
large scales, leading to deviations from the classical behaviour. 

Let me explain this point a bit further: In the underlying full LQG theory, the contributions 
to the (discrete) Hamiltonian operator depend on the state which describes the universe. As 
the volume grows (i.e., the universe expands), the number of contributions increases. Thus, the 
Hamiltonian constraint operator is expected to create new vertices of a lattice state (in addition 
to changing their edge labels), which in LQC results in a refinement of the discrete lattice |18j . 
Lattice refinement is also required from phenomenological reasons |1T|, I19j : for instance, it renders 
a successful inflationary era more natural [19J. The effect of lattice refinement has been modelled 
and the elimination of the instabilities in the continuum era has been explicitly shown. 

The appropriate lattice refinement model should be obtained from the full LQG theory. In 
principle, one should use the full Hamiltonian constraint and find the way that its action balances 
the creation of new vertices as the volume increases. Instead, phenomenological arguments have 
been used, where the choice of the lattice refinement model is constrained by the form of the 
matter Hamiltonian [20] . In particular, LQC can generically support inflation, and other matter 
fields, without the onset of large scale quantum gravity corrections, only for a particular model of 
lattice refinement [20J . This choice is the only one for which physical quantities are independent 
of the elementary cell adopted to regulate spatial integrations [21] . and moreover, it is exactly the 
choice required for the uniqueness of the factor ordering of the Wheeler-de Witt equation [22J. 



Allowing the length scale of the holonomies to vary dynamically, the form of the difference 
equation, describing the evolution, changes. Since the parameter hq determines the step-size 
of the difference equation, assuming the lattice size is growing, the step-size of the difference 
equation is not constant in the original triad variables. Let us consider the particular case of 

Hq ->■ fx (h) = fioV~ 1/2 , (20) 



suggested by certain intuitive heuristic approaches, such as noting that the minimum area used 
to regulate the holonomies should be a physical area [23], or that the discrete step size of 
the difference equation should always be of the order of the Planck volume. Moreover, this 
choice also results in a significant simplification of the difference equation, compared to more 
general lattice refinement schemes. The basic operators are given by replacing pq with /}. Upon 
quantisation [23] 

5^/2|/i) = er^l/i) ; (21) 
which is no longer a simple shift operator since jl is a function of p. Changing the basis to 

J H 3 

one gets 

e~ ifl ^\v) = e - m, ^\v) = \v + no) ■ (23) 

Thus, in the new variables the holonomies act as simple shift operators, with parameter length 
Hq. In this sense, the basis \v) is a much more natural choice than \p). One can then proceed 
as in the previous case of a fixed spatial lattice and write down the Hamiltonian constraint. 



(22) 



3.1. Constraints on inflation 

Let me consider first a fixed and then a dynamically varying lattice and solve the second order 
difference equation in the continuum limit. Two constraints can be imposed on the inflaton 
potential: the first one so that the continuum approximation is valid, and the second one so 
that there is agreement with the CMB measurements on large angular scales. A combination 
of the two constraints in the context of a particular inflationary model, will give us [20] the 
conditions for natural and successful inflation within LQC. 

More precisely, let us separate the wave-function ^(p, <j)) into ^(p, 4>) = T(p)<&((j)) and 
approximate the dynamics of the inflaton field, (f>, by setting V((f>) = V^p 5-3 / 2 , where is 
a constant and 5 = 3/2 in the case of slow- roll, to get [20J 
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with solutions [20] 



T (p) w p 



-(9+2<5)/8 



25 + 3 



7T 



C\ cos ( x 



3tt 



7T 



2(25 + 3) 4 



C<i sin ( x 



3tt 



7T 



2(2(5 + 3) 4 



_ (25) 
where x = 4^/5^(25 + 3)~V M+3)/4 and (3 = 96/(k^ 2 ). 

Without lattice refinement, the discrete nature of the underlying lattice would eventually 
be unable to support the oscillations and the assumption of pre-classicality will break down, 
implying that the discrete nature of space-time becomes significant on large scales. For the end 
of inflation to be describable using classical General Relativity, it must end before a scale, at 



which the assumption of pre-classicality breaks down and the semi-classical description is no 
longer valid, is reached. Let me quantify this constraint: The separation between two successive 
zeros of T (p) is 



Ap = -^P^ 25)/A ■ (26) 

For the continuum limit to be valid, the wave-function must vary slowly on scales of the order 
of fJL c = 4/2, leading to [20] 

Ap>^o(^f) 3/2 p~ 1/2 , (27) 



which implies [20] 



27vr 2 



v * < ^m^ m ■ < 28 > 

Assuming slow-roll inflation, we set 5 ~ 3/2. Setting /xo = 3\/3/2, the constraint on the 
inflationary potential, in units of h = 1, reads 

V(((>) & 2.35 x 10~ 2 / P1 4 , (29) 

which is a weaker constraint than the one imposed for fixed lattices, namely [20J 

V* « lO- 28 /^ 4 , (30) 

assuming that half of the inflationary era takes place during the classical era. 

One can further constrain the inflationary potential so that the fractional over-density in 
Fourier space and at horizon crossing is consistent with the COBE-DMR measurements, namely 

« 5.2 x 10-^ . (31) 

To do so, one must however adopt a particular inflationary model. As such, let us select 
V(cj>) = m 2 (j) 2 /2. Combining the two constraints we obtain [20] 

m <, 70(e" 2JVcl ) Mpi and m <, 10 M P i , (32) 



for the fixed and varying lattices, respectively. 

In conclusion, the requirement that a significant proportion of a successful (in particular 
with respect to the CMB measurements) inflationary regime takes place during the classical era, 
imposes a strong constraint on the inflaton mass. Such a constraint turns out however to be 
much softer, and therefore rather natural, once lattice refinement is considered. In this sense, 
we argue [20J that lattice refinement is essential to achieve a successful inflationary era, provided 
inflation can be generically (i.e., with generic initial conditions) realised. 



3.2. Constraints on the matter Hamiltonian 

Given a lattice refinement model, we will show [19] that only certain types of matter can be 
allowed. To do so, we parametrise the lattice refinement by A and the matter Hamiltonian 
by 5 and solve the Hamiltonian constraint. The restrictions on the two-dimensional parameter 
space will become apparent once physical restrictions to the solutions of the wave-functions are 
imposed [T9] . 

To be more precise, consider 

/2 = fi fx A , (33) 
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Figure 1. The allowed types of matter content are significantly restricted. For a varying 
lattice (A ^ 0) it is not always possible to treat the large scale behaviour of the wave-functions 
perturbatively (dashed line with crosses) ' 1T)|. 



leading to 

v = — n TV • (34) 

Being only interested in the large scale limit, we approximate |19j the matter Hamiltonian by 
Htj, = i> 5 e(4>), implying 

e = e = -v- s H g y . (35) 

A necessary condition so that the wave-functions are physical, is that the finite norm of the 
physical wave-functions, defined by L_^ dv\v\ s ^i^2, is independent of the choice of <j> = <fio. 
The solutions of the constraint are renormalisable provided they decay, on large scales, faster 

To solve the constraint equation, we need to specify the from of H^, which has in general 
two terms with different scale dependence. Being interested in the large scale limit, there is one 
dominant term, allowing to write [TH] 

m<p = e»(<f>)v S " , (36) 

where the function e v is constant with respect to v. 

Solving [19] the constraint equation, we only consider the physical solutions. The large 
scale behaviour of the wave-functions must be normalisable, which is a necessary condition for 
having physical wave- functions, while the wave- functions should preserve pre-classicality at large 
scales, which is a necessary condition for the validity of the continuum limit. We thus obtain [19J 
constraints to the two-dimensional parameter space (A, 8), shown in Fig. 1 for A in the range 
< A < —1/2, imposed from full LQG theory considerations |17j . 

In conclusion, the continuum limit of the Hamiltonian constraint equation is sensitive to the 
choice of the lattice refinement model and only a limited range of matter components can be 
supported within a particular choice. 



3.3. Uniqueness of the factor ordering in the Wheeler-de Witt equation 

I will now show that the lattice refinement model fi = fio^ 1 ^ 2 , argued [22] to be the only 
one achieved by both physical considerations of large scale physics and consistency of the 
quantisation structure, is also the only model which makes the factor ordering ambiguities of 
LQC to disappear in the continuum limit [22 ■. 

Let me be more specific. Writing the gravitational part of the Hamiltonian constraint in 
terms of the triad and the holonomies of the connection, one realises that there are many ways 
of doing so. Considering for example, 
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(37) 



there are many possible choices of factor ordering that could have been made at this point; 
classically, the actions of the holonomies commute. However, each of these factor ordering 
choices leads to a different factor ordering of the Wheeler-de Witt equation in the continuum 
limit. The action of the factor ordering chosen in Eq. (|37|) leads [22] to 



e ijk tr [hji^ 1 h^ x h k 
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(38) 



while other choices have certainly a different action. 

Defining V\v) = V v \v), with V the volume operator with eigenvalues V u , the action of the 
above factor ordering on a general state |^) = ^2 u ipu\ l/ ) in the Hilbert space, reads [22J 



e ijk tr (hihjhi 1 h- ] 1 h k 
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(39) 



the action of any other factor ordering choice can be obtained in a similar manner. By noting 
that the volume is given by 

v v \ v )~\n(y)\ z,2 \v) , (40) 

where n{y) is obtained by 



l-A 



»o(l-A) 



we find 



Kinfc 



(y db nk) a 



3/[2(l-A)] 



(41) 



(42) 



where a = no (1 — A) /k. 

We then take the continuum limit of these expressions by expanding ip u ~ ip (y) as a 
Taylor expansion in small k/v. For the particular factor ordering chosen above, the large scale 
continuum limit of the Hamiltonian constraint reads 1 22 1 : 



lim eyktr ( h\h\h- h- h k 
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One can easily confirm that one obtains the same continuum limit for the Wheeler-de Witt 
equation, only for the choice A = —1/2 [22], in which case the Wheeler-de Witt equation 
reads [22] 



lim C, 




Thus, there is only one lattice refinement model, namely jl = /io/i -1 / 2 , with a non-ambiguous 
continuum limit. 

In conclusion, phenomenological and consistency requirements lead to a particular lattice 
refinement model, implying that LQC predicts a unique factor ordering of the Wheeler-de Witt 
equation in its continuum limit. Alternatively, demanding that factor ordering ambiguities 
disappear in LQC at the level of Wheeler-de Witt equation leads to a unique choice for the 
lattice refinement model. 



3.4- Numerical techniques in solving the Hamiltonian constraint 

Lattice refinement leads to new dynamical difference equations, which being of a non-uniform 
step-size, imply technical complications. More precisely, the information needed to calculate the 
wave-function at a given lattice point is not provided by previous iterations. This becomes clear 
in the case of two-dimensional wave-functions, as for instance in the study of Bianchi models or 
black hole interiors. I will present below a method |24j based on Taylor expansions that can be 
used to perform the necessary interpolations with a well-defined and predictable accuracy. 

For a one-dimensional difference equation defined on a varying lattice, the Hamiltonian 
constraint can be mapped onto a fixed lattice simply by a change of basis |24j . This method is 
however not useful for the two-dimensional case, where the Hamiltonian constraint is a difference 
equation on a varying lattice [T7] : 
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where 6^ and 6 T denote the step-sizes along the \x and r directions, respectively. 

In the case of lattice refinement, <5 M and 6 T are decreasing functions of [i and r, respectively, 
and the data needed to calculate the value of the wave-function at a particular lattice site are 
not given by previous iterations, as it is illustrated in FigEl We propose |24] to use Taylor 
expansions to calculate the necessary data points. Let us assume that the matter Hamiltonian 
acts diagonally on the basis states of the wave-function, namely 
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Given a function evaluated at three (non-collinear) coordinates, the Taylor approximation to 
the value at a fourth position is 
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Figure 2. (a) For the fixed lattice case the two-dimensional wave-function can be calculated, 
given suitable initial conditions (solid circles), (b) In the case of a refining lattice, the data 
needed to calculate the value of the wave-function at a particular lattice site (open square) is 
not given by previous iterations (solid squares) 



where the Taylor expansion is taken about the position (x2,y2), we have defined Sfj = X{ — Xj 
and b\- = in — i/j , and the differentials can be approximated using 

f(x uyi ) = f^jrf + Sfr^ +% 2 ^-\ _.+••• , (50) 
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where the dots indicate higher order terms. 

For slowly varying wave- functions, linear approximation is very accurate and higher order 
terms in Taylor expansion can only improve the accuracy by 10 -2 % [23]. This method can be 
applied in any lattice refinement model, while its accuracy can be estimated. 

By using this Taylor expansions method, we confirmed [23] numerically the stability 
criterion of the Schwarzchild interior, found [i7\ earlier using a von Neumann analysis, and 
investigated |24] the way that lattice refinement modifies the stability properties of the system. 

Let me now show how the underlying discreteness of space-time leads to a twist |24j in the 
wave- functions, for both a fixed and a varying lattice model. In the case of a constant lattice, 
an initially centred Gaussian will move to larger fj,, as r is increased [24 . However, for regions 
where the lattice discreteness is important, this no longer holds, and the value at one lattice 
point introduces a non-zero component to the value at a lattice point with larger \i coordinate, 
i.e., the wave-function moves to larger \x. This implies the existence of some induced rotation 
on the wave-function due to the underlying discreteness of the space-time. Including lattice 
refinement, this effect persists [23]- I n other words, there is no motion for r>(5 T , while in the 
case of lattice refinement this requirement is reached for lower r, since 6 T reduces as r increases. 
As the wave-function moves into a region in which the discreteness of the lattice is important, 
a motion will be induced, as in the constant lattice case. 



4. Anisotropic LQC 

Various aspects of anisotropic cosmologies have been studied within LQC in the past [TTl [25]. 
however the first full and consistent quantisation of a Bianchi I cosmology (the simplest of 
anisotropic cosmological models) was achieved recently in Ref. [26]. Moreover, the link back 
to the underlying full LQG theory has been strengthened, by considering the flux of the triads 
through surfaces consistent with the Bianchi I anisotropic case. With the quantisation of the 
Bianchi I model under control, it is possible to ask whether LQC features, obtained within the 
context of isotropic Friedmann-Lemaitre-Roberston- Walker cosmology, are robust, at least with 
respect to this limited extension of the symmetries of the system. Bianchi type I models, apart 
their simplicity, they present a particular interest for space-like singularities within the full LQG 
theory. Following the same vein as for the isotropic case, a massless scalar field plays the role of 
an internal time parameter. In the absence of such a field, physical evolution can be found by 
constructing families of unitarily related partial observables, parametrised by geometry degrees 
of freedom [27] . 

4-1- Unstable Bianchi I LQC 

Studying stability conditions of the full Hamiltonian constraint equation describing the quantum 
dynamics of the diagonal Bianchi I model, I will show [28] that there is robust evidence of an 
instability in the explicit implementation of the difference equation. On the one hand, such a 
result may question the choice of the quantisation approach, the model of lattice refinement, 
and/or the role of ambiguity parameters. On the other hand, one may argue that such an 
instability may not be necessarily a problem since it might be that unstable trajectories are 
explicitly removed by the physical inner product. 

Consider the difference equation arising from the loop quantisation of the Bianchi I model [26] : 




ttG 



{v + 2) yjv + 4*+ (Ai, A 2 , v\ T) - {v + 2) y/V*+ (Ai, A 2 , u; T) 



(v - 2) (Ai, A 2 , v- T) + (v - 2) y/\u - (Ai, A 2 , u; T) ,(52) 



where 



(X 1 ,X 2 ,u;T) 



^2 *(oiAi,OjA 2 ,i / + 4;T) 

#j=(Q,l,2) 



^4 (Ai,A 2 ,i/;T) 



^(a i X 1 ,a j X 2 ,v-4:;T) 

^ = (-3-2,0) 



(Ai,A 2 ,^;T) 



* {aiXi,ajX 2 ,iy;T) 

iyy=(-l,0,l) 





(53) 




(54) 



Note that Ai, A 2 , A3 are related to the volume of the fiducial cell V by 




(55) 





Figure 3. The geometry of the points used in the difference equation that results from the 
Hamiltonian constraint, for the Bianchi I model [28J. 



with 7 = sgn(pip2P3)|7| and 

v = 2AiA 2 A 3 , (56) 

We will investigate the stability of the vacuum solutions, in which case the solution is static, 
namely ^ (Ai, A2, v,T) = (Ax, A2, v), simplifying considerably Eq. ([52]) . whose geometry is 
drawn in Fig. [3l To search for growing mode solutions we will perform a von Neumann analysis. 
In addition to specifying the boundary conditions on the v and u — 4 planes, we are also required 
to specify the value at five of the points given in ^\ (Ai,A2,^). There are in total 23 values 
that are required and with such initial data the difference equation, Eq. ()52j) . can be used to 
evaluate the 24 th point. Once this point has been evaluated, it can be used to move the central 
point and evaluate the wave-function at subsequent positions in the v + 4 plane. 

Following the standard von Neumann stability analysis, we decompose the solutions of the 
difference equation into Fourier modes and look for growing modes. Considering the ansatz [28J 

f (Ai, A 2 , v) = T (Ai) exp (i (wA 2 + xy)) , (57) 



where we have chosen the Ai direction to be the direction in which the v + 4 plane is evolved, 



the difference equation becomes 



e 4 X i £ T(a i Ai)e i(aJ ^ A2+ ^ ) 

i#=(0,l,2) 



+ 



^•=(-1,0,1) 
i/- 2' 



i/ + 4 



i/ + 2 / V ^ + 4 



u-2\ \v-A 



I/ + 2/V ^ + 4 



^ T(a,Ai)e i( ^ A2+Xiy) 

.^•=(-2,0,3) 

4x * T(a J Ai)e i( ^ A2+ ^ ) 
^=(-3,-2,0) 



(58) 



As it has been explicitly shown in Ref. [28J, expanding in terms of small 1/v, the difference 
equation can be written in the form of a vector equation as 



MiT 3 = M 2 T 2 



(59) 



where we have defined the vectors 
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r(oj_iAi) 

T(a i - 2 \ 1 ) 

r(oi_ 3 Ai) 

T(ai_ 4 Ai) 
T(ai_ 5 Ai) 
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2,3 



(60) 



and the matrices 
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The condition for stability can be written as follows: If 
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(62) 



where A are the eigenvalues of the matrix (Ml) - M 2 , the amplitude T (a^Xi) is less than that 
of previous points, in other words, the difference equation is stable. Let us define the parameter 
A = io\ 2 jv. The presence of an instability in the difference equation has been demonstrated [28 
in several ways: (i) there is a particular set of critical modes, A = (2n — l)ir/2, with n G Z, 
for which the system is unstable; (ii) in the large v limit, the system is unstable for the modes 
A = 7r/4 and x = 0; (iii) the system is unstable for a general u, for modes that approach the 
critical value. 

In conclusion, the difference equation, Eq. (j52[) . is unconditionally unstable, in the sense that 
there is no region of (Ai, A2, v) in which Eq. (|52|) is stable. 



4-2. Lattice refinement from isotropic embedding of anisotropic cosmology 

Given a consistent quantum anisotropic model, one can find |29| isotropic states for which the 
discrete step-size of the isotropically embedded Hamiltonian constraint is not necessarily that 
of the p~ l l 2 (new quantisation) step-size. The choice of different embeddings has important 
consequences for the precise form of discretisation in the isotropic sub-system. In this sense, 
lattice refinement could be interpreted as being due to the degrees of freedom that are absent 
in the isotropic model [29 . 

In the standard approach, isotropic states are taken to be those in which the three scale 
factors are Ai = A2 = A3, or more precisely, defining the volume of the state as v = 2A1A2A3 to 
eliminate one of the directions (A3), the map 

I* (Ai, A 2 , „)) | J] * (Ai, A 2 , v)) = |tf (*/)) , (63) 

Ai,A2 

produces isotropic states. Working with the three scale factors A,, one can show |29j that there 
is an ambiguity in exactly what the volume of such an isotropic state is. Consider the state 

I* (Ax, A 2 , A 3 )) = j (Ai, A 2 , A 3 )) + |tf (A 3 , Ai, A 2 )) + I* (A 2 , A 3 , Ai))] , (64) 
with the restriction 

(* (Ai, A 2 , A 3 ) |* (A 3 , Ai, A 2 )) = (tf (Ai, A 2 , A 3 ) |tf (A 2 , A 3 , Ai)) = (* (A3, A x , A 2 ) |tf (A 2 , A3, A x )) , 

on the anisotropic states. The expectation values of the scale factors along each direction of 
such a state are 

/\ \ Ai + A 2 + A 3 

(^) = g • (65) 

The measured scale factor is equal in each direction and is given by the average of the scale 
factors of the underlying, anisotropic states. However, the measured volume of such a state is just 
v = 2AiA 2 A3, which is not the cube of the measured scale factor. Thus, while it is the eigenvalue 
of the anisotropically defined volume operator, it is not necessarily what we would measure as the 
volume. Essentially, the reason for this is that while the scale factors Aj are measured to be equal 
in each direction, they are not eigenvalues of the state, i.e., Xi\^ (Ai, A 2 , A3)) 7^ X{\^ (Ai, A 2 , A3)). 
However, both the average and the product (i.e., the volume) of the scale factors are eigenvalues. 
It is this ambiguity, that leads to the possibility of deviations from the standard isotropic case. 

In conclusion, the difference between the two procedures is essentially due to what one 
considers to be more fundamental, the volume of the underlying states (v) or the measured 
volume of the symmetric state ((A) 3 ), which are not necessarily equal. Choosing v leads to the 
new quantised Hamiltonian of isotropic cosmology, while choosing (A) 3 results in some kind of 
different lattice refinement. This lattice refinement is significantly more complicated that the 
single power law behaviour, yr 4 , usually considered. 



5. Conclusions 

Loop Quantum Gravity proposes a method of quantising gravity in a background independent, 
non-perturbative way. Quantum gravity is essential when curvature becomes large, as for 
example in the early stages of the evolution of the universe. Applying LQG in a cosmological 
context leads to Loop Quantum Cosmology which is a symmetry reduction of the infinite 
dimensional phase space of the full theory, allowing us to study certain aspects of the theory 
analytically. The discreteness of spatial geometry, a key element of the full theory, leads to 
successes in LQC which do not hold in the Wheeler-de Witt quantum cosmology. 

By studying phenomenological consequences of LQC we can get some useful insight about the 
full LQG theory and get an answer to some long-standing questions in early universe cosmology. 
Here, I have briefly described some of the phenomenological aspects of LQC. 
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